Period-doubling bifurcation in strongly anisotropic Bianchi I quantum cosmology 



Michael Bachmann* and Hans-Jiirgen Schmidt^ 
Institut fur Mathematik, Universitat Potsdam, PF 601553, D-14415 Potsdam, Germany, 
Institut fur Theoretische Physik, Freie Universitat Berlin, Arnimallee 14, D-14195 Berlin, Germany 

(February 7, 2008) 

We solve the Wheeler-DeWitt equation for the minisuperspace of a cosmological model of Bianchi 
type I with a minimally coupled massive scalar field (j> as source by generalizing the calculation of 
Lukash and Schmidt [1]. Contrarily to other approaches we allow strong anisotropy. Combin- 
ing analytical and numerical methods, we apply an adiabatic approximation for (f>, and as new 
feature we find a period-doubling bifurcation. This bifurcation takes place near the cosmologi- 
cal quantum boundary, i.e., the boundary of the quasiclassical region with oscillating t/j- function 
where the WKB-approximation is good. The numerical calculations suggest that such a notion of 
a "cosmological quantum boundary" is well-defined, because sharply beyond that boundary, the 
WKB-approximation is no more applicable at all. This result confirms the adequateness of the 
introduction of a cosmological quantum boundary in quantum cosmology. 
PACS 98.80; 98.70 

I. INTRODUCTION 

The idea to consider the whole Universe as one quantum system has already attracted many researchers. One of 
the approaches to give this idea a physical meaning is to consider the superspace S of all possible spatial 3-geometries 
together with the matter field degrees of freedom as supercoordinates, and the supermomentum space S* is formed 
from the corresponding second fundamental forms together with the momenta of the matter fields. The dynamics 
will then be defined by the Einstein field equation accompanied by the matter held equations. In the last step, this 
system will be quantized. 

Of course, one has to deal with an infinite number of degrees of freedom already from the geometric part. The, 
up to now only, way to tackle such a system is to linearize in almost all degrees of freedom and to take the non- 
linearities only from the remaining finitely many ones. But even this system is almost untractable. So, the idea of 
the minisuperspace arose: One restricts to a finite-dimensional space of spatial geometries and to matter fields with 
finitely many components, disregarding or even ignoring all other degrees of freedom. 

From the first glance one could believe that such a big simplification would lead to a picture which has nothing to 
do with the real Universe's evolution. However, the results given in the last years are encouraging: The dynamics 
of the minisuperspace can be written as an equivalent mechanical system, and then the Schrodinger equation for 
this system carries the name Wheeler-DeWitt equation which can be solved for simple systems and give already a 
surprisingly good picture of the evolution, even for the case that the geometric part is restricted to only one degree 
of freedom, the "radius" of the Universe. This corresponds to the isotropic Friedmann Universe model. For this case, 
closed-form solutions for the Wheeler-DeWitt equation exist. 

In the ordering with respect to simplicity, the next possible geometry is the Bianchi type I model, where the spatial 
inner geometry is flat, but the expansion is allowed to be anisotropic. For this case, we have 3 degrees of freedom 
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for the geometry: the expansion rates into the 3 spatial directions.^] For this case, the Wheeler-DeWitt equation is 
already quite complicated, so one uses an approximation to solve it. One of the most powerful of these approximation 
schemes is the Wentzel-Kramers-Brillouin (WKB)-approximation. In Rcf. [1], the first order WKB-approximation for 
the massive scalar field in the Bianchi type I model has been deduced, but - as far as we are aware - the higher order 
WKB-approximations have not yet been calculated up to now for this model. The analogous calculations for Bianchi 
type IX have been done by Amsterdamski [2]. In both cases, the anisotropy degrees of freedom had been assumed 
to be small. In the present paper we allow also large anisotropies. Ref. [3] deals with a similar model allowing large 
anisotropics, too; however, the authors of [3] use a simplicial minisuperspace model, so the results are not directly 
comparable to ours. 

The text is organized as follows: Sect. II presents the geometry of the Bianchi type I cosmological model, Sect. Ill 
the corresponding Whccler-DeWitt equation and Sect. IV its solutions both analytically and numerically. Sect. V 
shortly reviews Refs. [6 - 42], i.e., the earlier developments of the topic, and discusses the results. 

In the appendix, we show how the different versions to solve the factor-ordering problem influence the solutions of 
the Wheeler-DeWitt equation. 



II. THE BIANCHI TYPE I MODEL 

The metric of a spatially flat cosmological model is deduced as follows: One requires that an abelian 3-dimensional 
isometry group acts transitively on 3-dimensional spatial hypersurfaces which consequently have to be flat 3-spaces. 
This restricts the possible metrics to 

ds 2 = N 2 dt 2 - g ij dx i dx j (2.1) 

where N — N(t) is the lapse function which can be put to N = 1 by a time reparametrization, i,j — 1, 2, 3, and gij 
represents a symmetric positive definit matrix whose components depend on t only. So we have 6 free components 
from the first consideration. However, restricting to metrics being a solution of Einstein's field equation we can 
simplify as follows: At any initial time, say t = 0, we can choose the initial condition 

g tj = 5ij = diag(l, 1,1) at t = (2.2) 



without loss of generality, and then by a spatial rotation not changing Eq. (^2), the second fundamental form, 
i.e. dgij/dt, can be brought into diagonal form. By use of the Einstein equation one can show that under these 
circumstances, g^ will keep its diagonal form for all times. Thus: The Bianchi type I model containing non-diagonal 
terms does not represent a generalization, and we can say: Without loss of generality let 



g i3 = di ag (A 2 (t),B 2 (t),C 2 (t)) (2.3) 



with certain positive functions A, B, C. 



After an obvious rearrangement of the terms we get now from Eqs. (2.1) and (2.3) 

ds 2 = N 2 dt 2 - e 2a (e 2s+2 ^ r dx 2 + e 2s - 2 ^ r dy 2 + e^dz 2 ^) (2.4) 



1 Strictly speaking, the axially symmetric Bianchi type I model, where 2 of the 3 scale factors coincide, is even simpler, but as 
one will see from the calculations: restricting to axial symmetry does not really simplify the procedure. 
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where a, r, s represent 3 arbitrary real functions depending on t. The Hubble parameter is the mean expansion, i.e. 
with a dot denoting d/dt we get 




. In the notation of Eq. (2.4), H — da/dt is the Hubble parameter, and we restrict our considerations to the expanding 
Universe, i.e., to H > 0. The model goes over to the spatially flat isotropic Friedmann model if the dimensionless 
anisotropy parameter 



l-Vr 2 +s 2 (2.5) 



11 H 
tends to zero. 

It is essential to observe that this approach leading to the Wheeler-DeWitt equation breaks the Lorentz invariance 
of the system already at this level of geometry: The 3+1 decomposition of space-time is already made from the 
beginning. The remaining degrees of freedom within the geometry are the following: t(t) as arbitrary function as 
long as N{t) remains unspecified, and t(t) = ±t after having fixed N = 1. Permutations of the 3 spatial coordinates, 
and for the axially symmetric case also spatial rotations: If r = 0, then a rotation in the xy-plane is an additional 
symmetry. We presented this geometric part so explicitly, because contradicting statements exist about this behavior 
in the literature; example: In our interpretation, ds 2 = dt 2 — dx 2 — dy 2 — dz 2 and ds 2 = dt 2 — t 2 dx 2 — dy 2 — dz 2 
represent different geometries inspite of the fact that they are locally isometric space-times. 

We use a one-component real massive scalar field 4> as source, and we interpret it as follows: Either it can really 
be such a scalar field (i.e. a spin zero field), e.g. a Higgs field which dominated the early Universe's evolution but 
disappeared after a symmetry breaking effect. Or, it mimics any realistic matter field in a region where spin is 
negligible. Let m denote that mass of the scalar field. For the GUT-mediated inflationary model, one considers the 
following order of magnitude: m w 10~ 5 mpi, where mpi is the Planck mass.0 

Up to now, we did not restrict the spatial coordinates x, y and z. At the classical level (i.e. for locally solving the 
Einstein field equation), it makes no difference, whether they cover all the reals or whether they are cyclic ones. But 
already for global classical considerations it makes a difference: Supposed, x, y and z are cyclic, then even for the 
case r = s = the metric is not spatially isotropic, because the natural length of a closed spatial geodesic depends 
on the chosen direction. However, we get the essential difference only after quantization. If nothing different is said, 
we assume x, y and z to be cyclic with modulus 1, i.e. x and x + 1 represent the same point etc. 



III. THE WHEELER-DEWITT EQUATION 

The Wheeler-DeWitt equation is always a zero energy Schrodinger equation due to time-reparametrization in- 
variance of the gravitational action. However, this is not really a problem, because the Schrodinger equation with 
non-vanishing energy can be brought to zero-energy form by a suitable shift of the potential, see e.g. Ref. [4] for a 
class of transformations of this type. 

In the present paper, we do not consider the entropy of the system; we only want to mention the recently published 
result of Kleinert [5] on how the incorporation of entropy can change the model. 



2 In units with c = 1, the Planck mass is about 10 19 GeV. 
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To formulate the equations that lead to the Wheeler-DeWitt equation [6] as easily as possible, it proves useful to 
apply units such that Ti = c = AirG/S — 1. Then the Lagrangian of the minimally coupled scalar field in Einstein's 
theory of gravity (see the first of Refs. [1] for details) is given by 



L = L g + L m , (3.1) 



where the gravitational part is 



and the matter Lagrangian is 



L g = ^{r 2 +S 2 ~a 2 ) (3.2) 



Lm - - 



^e 3 «(0 2 -mViV 2 ). (3.3) 



Here we supposed the lapse function N to be constant in time. The minisuperspace is 4-dimensional, its coordinates 
are = (q°, q 1 , q 2 ,q 3 ) with q° — a, q l — <j>, q 2 = r, and q 3 = s. The usefulness of the Misner parametrization 



Eq. (2.4) of the anisotropy is obvious: The kinetic part of the Lagrangian is already in diagonal form. By defining 
the minisuperspace metric = e 3a diag(— 1, 1,1,1) the Lagrangian (3.1) is compactly rewritten as 



2N 

Thus the introduction of canonical momenta 



^(/,^-mV) 2 ^ 2 e 3 *°). (3.4) 



P " = W = N f ^ q (3 - 5) 
allows us to write down the classical Hamiltonian of the minisuperspace 

U = f + ™V) 2 e 39 °) = °- ( 3 - 6 ) 

Time-reparametrization invariance implies ^ — 0, i.e., Ti — 0. The quantization in coordinate space is done as usual 
by going over to the appropriate operators p M = —ihd/dq^ with the partial derivative d. 

From now on, we will choose two different versions for the operator ordering in p^p^: 

First variant, variant A, which is the most simple one, follows the interpretation in Ref. [1] that the solutions of 
the forthcoming Wheeler-DeWitt equation will not essentially depend on this choice. 

Second, variant B, solves the factor ordering problem by applying the supercovariance principle: It says, that the 
essential equations have to be covariant ones with respect to the (mini-)superspace metric. Taking variant B, we have 
to replace p M by p^ = —ihV /V q^ where V denotes the supercovariant derivative^ defined by /^„. 

To cover both variants simultaneously, we introduce a parameter e as follows: e = for variant A and e = 1 for 
variant B. Accordingly, we denote the two variants of the D'Alembert operator^ O e by 

D ° = -^^- ^ 

and 



3 Here we apply the fact that the metric is always covariantly constant, so different orderings give the same result. 
4 The minus sign in eqs. (3.7/3.8) is inserted to compensate for the factor i 2 in front of p /J p /J '. 
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For the world function ipiq^) we get the relations 
where^J 

a tjj = e~ 3a [ip^o - - ^,22 - ip, 



(3.8) 



(3.9) 



(3.10) 



The stationary Schrodinger equation with zero energy in minisuperspace TLi]j{q^) = with the wave function ^(q^) 
and the Hamilton operator TL is called Wheeler- DeWitt equation and possesses the form 

N 



Win = - (h 2 U t + mV) 2 e 3< ?°) ^) 







for the Bianchi type I Universe (2.4). In details, we get after dividing by Ne /2 







<9 2 



d d 2 d 2 d 



+ mVe ba ip(a,<f>,r,s). 



(j) 2 dr 2 ds 2 

In what follows we derive an approximation scheme to solve the Wheeler-DeWitt equation. 



(3.11) 



(3.12) 



IV. SOLUTIONS OF THE WHEELER-DEWITT EQUATION 

Before we go over to solve Eq. (3.12) we want to clarify in subsection A what happens if we neglect the anisotropy 
from the beginning. In subsections B to D we analytically obtain solutions in a well-defined approximation, and in 
subsection E we present the results of our numerical calculation. 

A. Isotropic models 



For variant A, i.e., e = in Eq. (3.11) we get 



= ft 



d 2 



d 2 



da 2 dcf) 2 



+ m 2 <j> 2 e 6a ) i>{a,(t>) 



(4.1) 



i.e., we simply remove all r- and s-dependencies. For variant B, however, this recipe does not work because the 
covariant derivatives mix the coordinates. After a short calculation one finds out that here also Eq. (4.1) is the 
correct Wheeler-DeWitt equation. The reason for this coincidence is as follows: The 2-dimensional superspace is 
conformally flat, and the DAlembertian is conformally invariant, cf. the second paper in Ref. 1 for more details. 



Defining a — e a > as new variable we have to replace by cl-§^, and we get from Eq. (4.1) now 

d 2 



0= ft 



d_ d_ 

da da dcf) 2 



+ m 2 (/) 2 a b ip(a, </>) 



(4.2) 



a version of the Wheeler-DeWitt equation favoured in Ref. [7] 



5 If the dimension of the minisuperspace is D instead of 4, then the factor 3 in front of V'.oin Eq. (3.9) will be replaced by 
3(D - 2)/2, so D and Di coincide for D = 2. 
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B. Separation ansatz 



Now we look for the solutions of Eq. (3.12). First, we make a separation ansatz 

tp(a, <f>, r, s) = X (a, <j>) ■ p{r) ■ a(s). (4.3) 
It turns out that ppij p = ci and ct 33 /cr = c 3 represent constants, i.e., we get 

p{r) = pi cxp( A /c 2 ~ r) + p 2 cxp(-^/c2 r) 

and 

a(s) = ci exp(^/c^ s) + cr 2 cxp(-^ci" s) 

with arbitrary constants pi and Oi. For negative values C2, the function p(r) has a sinus-function behaviour. 

One can interpret these cases as follows: If both c 2 and C3 are negative, then we have plane waves in r- and 
s-direction. This is in agreement with the fact, that a translation into r- or s-direction can be compensated by a 
coordinate transformation in metric (2.4) (i.e., by multiplying x, y, and z with suitable constants), so that all r- and 
s-values should be equally probable. However, this interpretation works only in case that we allow x, y, and z to 
cover all the reals. If we, however, restrict x, y, and z to be cyclic coordinates, then this argument in favour of equal 
distribution of the r- and s- values is no more valid. 

If, on the contrary, one of the constants c 2 or C3 is non-negative, then (up to singular exceptions) the product 
p(r) ■ a(s) tends to ±00 as r or s docs. This allows the following interpretation: The probability to have small 
anisotropy is exponentially small. 

For the moment we keep the further interpretation open and continue the calculation. We define the constant C\ 
via ci + c 2 + C3 = and insert Eq. (4.3) into (3.12); we get: 







d 2 
da 2 



3e 7i 

da 



d 2 



b 2 e 6a 



X(a,<p). 



(4.4) 



This equation is for variant A just the Whcclcr-DcWitt equation (4.1), however, now not the zero-energy equation 
but the equation with energy proportional to the constant c\ representing the anisotropic degrees of freedom. 

We redefine the coordinate a to v = a 3 = e 3a . The coordinate v is proportional to the spatial volume. Then Eq. 
(4.4) goes over to (after dividing by 9) 



d_ d_ d d 2 

dv dv dv d(j) 2 



(4.5) 



As one can see in comparison with Eq. (4.2): Instead of the a 6 term we now have a w 2 -potential. 



C. Adiabatic scalar field 

In this subsection, we assume the scalar field (f> to be almost constant. Then 7 defined by j 2 — m 2 (f> 2 /(9h 2 ) is an 
adiabatic constant. In Eq. (4.5) we omit the ^-derivative corresponding to the adiabatic approach. To this end we 
consider the following equation for the exact dependence on v of a new function \ 

+ Cl h 2 /9 + 7 V?i 2 ) X (v, <P) - E(<f>) ■ X (v, 4>), (4.6) 

where E{<j>) is an yet undetermined eigenvalue; here, <p plays the role of a parameter only. We define a new adiabatic 
constant A by A 2 = E((/))h^ 2 - ci/9. Then Eq. (4.6) becomes for e = 



d d d 
av ov ov 
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(4.7) 



Eq. (4.7) is a Bessel type ordinary differential equation. If we replace in Eq. (4.7) 71* by x and x(w, (j>) by y(x,4>) we 
get exactly Bessel's form 



r 2 

dx 2 



X ^L + (x 2 -A 2 ) y = 0, 
ax 



whose solutions are (the constants Cj depend on A) 

y = C 1 J a {x) + C 2 J-a{x) . 
The Ja are called cylinder or Bessel functions and have the development 



Mx) = f>i) fc (| ) A+2k (kiy 1 (r(A + k + 1))- 1 



k=Q 



Now we make the ansatzQ 



X{v,4>)= / b A ((/>)xA(v,<l>)dA, 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



where we assigned to \ the subscript A to emphasize xa being an eigenfunction with eigenvalue A. Thus the general 
solution of Eq. (4.7) is written as 



Xa(«, <f) = CiJa{jv) + C 2 J-a(iv) 
Inserting the ansatz (4.11) into (4.5) we get 

d 2 







-tr 



P>a(0)xM)] + &a(0 



lr j v—v— + — 

dv dv 9 



„2i2 
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x(v,(f>)dA. 



(4.12) 



(4.13) 



In the second term of Eq. (4.13) we replace the term in the brackets [ ] by the eigenvalue E((f>) according to Eq. (4.6). 
E(cj>) depends on C2, C3 and A only, i.e., there is no independence, and we write E instead of E(cf>) 







(4.14) 



In Eq. (4.14) we first write out the applied to the product, but we assume that only ^3-&a(</>) is essential, and 
disregarding the other terms^] and assuming that the x~ s are all independent we obtain the simple differential equation 



Ft 2 



b A {</>) = 



which possesses solutions 

b A (<t>) = B 1 e V ^ 4 ' /h + B 2 e- V Z' t > /h . 
Finally, we get the general wave function 



(4.15) 



(4.16) 



If only a discrete set of values A appear, Eq. (4.11) will be replaced by an analogous sum. 
This is usually called adiabatic approximation with respect to cf>. 
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V'cfL ("> & r > s ) = (Pi ex P {v^ r} + p 2 exp {- 



exp | A - (c 2 + c 3 
x [CiJ A (m0e 3Q /3?i) + C 2 J-a (m(f>e 3a /3fi)] dA 



c 2 r}) (cri exp {y^ s} + cr 2 exp {--y/c3 s}) 
)J9^+B 2 exp{-v/A-(c 2 +C3)/9 0} 



(4.17) 



with real and for the moment continuous eigenvalues c 2 ,C3. The superscript (4>) at ip indicates that we have treated 
the scalar field adiabatically. This solution has to be specified by appropriate choices for the constants pi, p 2 , cri, cr 2 , 
Bi, £> 2 , Ci, and C 2 from boundary conditions. Before discussing the result ( 4.17 ), however, we calculate the wave 
function by considering the scale factor a as the adiabatic variable. 



D. Adiabatic Scale Factor Approach 



Another adiabatic method to solve the Wheeler-DeWitt equation (4.4) is to treat the scale factor a as a slowly 
varying variable. 

Once more, we reexpress the scale factor a by v = e 3a , implying that we consider the Wheeler-DeWitt equation (4.5) 
in the following. Considering v as the adiabatic variable means that we neglect derivatives with respect to v in a first 
step. Then Eq. (4.5) becomes 



if 



1 



+ — ci + - m <j> v x(v, <fi) = E(v) x(v 



with E(v) being a still undetermined eigenvalue. By defining 



1 

lo = — mv, x 
3 



(4.18) 



(4.19) 



the differential equation ( 4.1S| ) is simply transformed into that for the dimensionless harmonic oscillator 

y(x) = 0. (4.20) 



By substituting x = x 2 and then transforming y(x) = e x / 2 y(x), we eventually obtain Kummer's differential equation 



d 2 y{x) 



+ (M ~ x) 



dy(x) 



vy{x) = 0, 



dx 2 ' ' dx 

where in our case fj, = 1/2 and v = (1 — rj)/&. This differential equation possesses the general solution 

y(x) = A x Fi{v, /it; x) + B x Fi{v - fi + 1, 2 - fi;x) x 1 ^, 



(4.21) 



(4.22) 

expressed with the help of the confluent hypergeometric function i-Fi(a, b;x). A and B are constants. Thus the partial 



problem of finding the solution of Eq. fl4.1q ) is done and the result is written as 

Xv (v^) = e-^ 2 / 2h A 



iFi ((1 - 77)74, 1/2; u(v)4> 2 /h) + B 1F1 ((3 - ry)/4, 3/2; uj(v)4, 2 /h) 



(4.23) 



with ^ denoting the eigenfunction to the real, continuous eigenvalue rj. 

For further processing to find an adiabatic solution for the Wheeler-DeWitt equation (4.5) we make the integral 
ansatz 



X{v,4>)= I g r ,{v)xr,{v,4>)drj. 

Jr) 



(4.24) 
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Inserting this into Eq. (4.5), we obtain (for e = 0) 



d d 

0= / h2v 7T v 7r [9ri(v)Xr,(v,<t>)} +9 v (v) 



_-I (/l \, +m W) 



X v (v,<j))dTi. 



(4.25) 



(4.26) 



dv dv 

The expression enclosed in the brackets [ ] in the second term is substituted by the eigenvalue 

E ( v ) = g [3ftmr]v - h 2 (c 2 + c 3 )] , 

following from Eq. (4. IS) and the definitions ( f4.19| ). In contrast to the preceding approach with an adiabatic scalar 
field, E(v) depends explicitly on the adiabatic variable which is in this case v. Note that E(v) also depends on 
the eigenvalues c 2 , C3, and 77. Now we utilize the assumption of adiabaticity with respect to v by neglecting terms 
containing derivatives dx v /dv and d 2 \Tj/9v 2 which appear in the first term of Eq. (4.25). Furthermore supposing 
that all functions Xv( v i 4>) are independent with respect to 77, each integrand in Eq. ( 4.25| ) vanishes. Thus we remain 
with solving the differential equation 



-2 d d 



hv o- v v d- v +Eiv) 



g(v) = 



or, more explicitly, 



2 



d 



v "o - 2 + V ~5T + KlV 

av z av 



g(v) = 



(4.27) 



(4.28) 



with K\ — mr)/3h and k\ = 4(c2 + c^)/9. Applying the transformation x = 2^/k\v and denoting the solution as 
y(x) — g(v(x)), Eq. ( 4.28 ) takes the same Bessel form as given in Eq. (4.8), whereas the index is K\ instead of A here. 

(4.29) 



For this reason the general solution of Eq. ( 4.28 ) reads 

g v (v) = Di J 4 (c 2 +c 3 )/9 (V 4m W 3?i ) + D 2 J-4(c 2 +c 3 )/9 (V 4m W 3?i ) 



with constants D\ and D 2 - 

Thus, the complete general wavefunction obtained with an adiabatic scale factor approach (indicated by super- 
script a) is found to be 

tp^l s (a,(t),r,s) = (piexp{^ci r} + p 2 exp{- % /c 2 " r}) (criexp{^ci' s}) + a 2 exp { - y 7 ^ s}) 

A 1 F 1 ((1 - r?)/4, 1/2; m#e 3a /3fi) + B J ' ^e 3a/2 1^1 ((3 - r,)/4, 3/2; m <b 2 e 3a /3h) 

with constants p\, p 2 , <Ti, <j 2 , A, B, D\, and D 2 to be determined via appropriate boundary conditions. The anisotropy 
quantum numbers c 2 and C3 are supposed to be real at this stage. Certain boundary conditions, such as requiring the 
wavefunction to vanish for infinitely large values of a and/or —a, lead to replacing the integral by a sum over r\. 



3)/9 ( v / 4mr ? e 3Q /3?i) + D 2 J_ 4(c2+C3)/9 U/Amr}e 3a /^n\ exp {-m^ 2 e 3Q /6fi} 

drj (4.30) 



E. Visualization of the world function 



Now we visualize the results of our numerical calculation of the world function (4.17). To get an impression about 
the shape of the probability amplitude, we concentrate on the absolute value of the wave function \ip\ which is given 
in these six diagrams in dependence of a, the logarithm of the cosmic scale factor, and the scalar field (p. Each 



figure is printed at fixed values of the other parameters in Eq. ( 4.17 ). We chose p\ = o\ = B\ = C 2 = and 
P2 = &2 = B 2 = C\ = 1 to ensure convergence for small and large values of a and (f>. The absolute wave function 
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\ip(ct, (j>) \ is not normalized, however, when comparing different figures (for instance Fig. [I] and Fig. g), higher function 
values mean an increased probability density. 




1.5 

FIG. 3. The same as Fig. [l], now ci = —4 and C3 = —7. The small discontinuities in the figure in the lower left corner 
may be numerical artefacts, so we do not comment them further. In comparison to Figs, [l] and [], the different spacing of the 
|?/;|-axis is still greater. 




FIG. 4. The wave function in the isotropic limit, represented by C2 = and C3 = 0. 
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3 

FIG. 5. The wave function with quantum numbers C2 = 1 and C3 = 1. This figure is also an example for the influence of 
the anisotropy upon the shape of the wave function. 




FIG. 6. Magnification of a subrange of Fig. [y. From the description it becomes clear which part of Fig. || it represents. We 
show it here to make clear, how quickly the WKB-approximation becomes a good approximation. 

For all these different ranges of parameters leading to the results shown in Figs. |l|-^ we find essentially the same 
behaviour: The right part of the figure is the region where the WKB-approximation is valid, that means, the solution 
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is close to a sinus-shaped wave with a slowly varying amplitude (see Fig. ^). The left hand side is the quantum region 
where no oscillations exist at all. The range where these two behaviours go into each other, is relatively sharply 
defined. This range is what one usually calls the "cosmological quantum boundary" , and these pictures show up that 
this notion is relatively well-defined. 

An interesting new feature of these figures is the following: They show a period-doubling bifurcation of the frequency 
if one looks from the bottom to the top of these figures with increasing values of (p. 

V. DISCUSSION 

Before we discuss the results of the present paper we give a short review to other papers: 

The original paper [6] by DeWitt was seminal to the whole development of quantum gravity. Its application to 
minisuperspace cosmological models has been developed by several researchers, e.g. the authors and editors of refs. 
[7] and [8]. At that time, the main discussion dealt with the closed isotropic Friedmann models, and the spatially flat 
Friedmann models had been included as less interesting limit cases, too. 

In [9] and [10], Schon and Hajicek perform a quantization of systems with quadratic constraints and discuss the 
Wheeler-DeWitt equation. In ref. [11], Landsman gives details which Hilbert space might be appropriate for the 
isotropic minisuperspace quantization, a topic, which was not much discussed before. 

In the recent preprint by Kim [12], the minimally coupled massless scalar field in an open isotropic Friedmann 
model has been discussed and its Wheeler-DeWitt equation be solved from the point of view that the Universe is 
created quantum mechanically from "nothing" . 

Recently, Capozziello and Lambiase [13] investigated the connection between the Hartle criterion for selecting 
correlated regions in the configuration space of dynamical variables and an associated Noether symmetry. This 
relationship serves to classify solutions of the Wheeler-DeWitt equations in semi-classical minisuperspace models. 
Thus, the oscillatory behaviour of a subset of solutions entails the presence of Noether symmetries which, in the 
consequence, select classical Universe models. 

The paper [14] by Schunck and Mielke is related to the models discussed here as follows: They apply the Wagoner- 
Bekenstein-Starobinsky-transformation for classifying inflationary solutions with scalar field as source, and this clas- 
sification should also apply to the corresponding Wheeler-DeWitt equations. 

In [15] and [16], a one-parameter set of minisuperspace metrics in arbitrary dimensions is considered, from which 
we have chosen only that one which gives classically the correct correspondence to Einstein's theory. In [15], the 
signature of the superspace metric in dependence on the signature of the underlying manifold is evaluated with the 
result that the normal-hyperbolic character of the Wheeler-DeWitt equation exists only for the Euclidean and the 
Lorentzian signature of the underlying manifold. 

In [17], Horiguchi, Maeda and Sakmaoto perform an expansion of solutions of the Wheeler-DeWitt equation in 
powers of the Planck length. Vilcnkin [18] compares several approaches to quantum cosmology, Kim and Page [19] 
discuss quantum Friedmann models and power-law inflation, Kim [20] compares quantum Friedmann models with 
conformally and minimally coupled scalar fields. In Ref. [21], Bleyer and Ivashchuk discuss multidimensional cosmo- 
logical models and their corresponding Wheeler-DeWitt equations. Page [22] solves the Wheeler-DeWitt equation for 
scalar fields as source. 

Ref. [23] represents the famous paper in which the "Hartle-Hawking boundary conditions" for the Wheeler-DeWitt 
equation have been derived. One takes a path integral over all such space-times V4 whose boundary is the prescribed 
spatial hypersurface V3. 
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In Ref. [24], Kiefer constructs wavepackets in minisuperspacc for a Friedmann Universe. An adiabatic approach is 
used in the case of a massive scalar field, thereby assuming the scale parameter a to be slightly changing only. In 
Zeh [24], these solutions are discussed under the point of view of the definition of the direction of time. Our Eq. 
(4.1) coincides with the form discussed in Ref. [24] (Eq. (6.5) of Zeh, and Eq. (2.2) of Kiefer). The approximation 
of our subsection IV D was used by Kiefer [24] for solving the Wheeler-DeWitt equation of a Friedmann Universe. 
In contrast to Kicfer's procedure, where the harmonic oscillator has discrete eigenvalues and thus the wave function 
decreases for large values of the scale factor a, we allow all real eigenvalues. 

Conradi [25] solves the Wheeler-DeWitt equation for Bianchi type IX and a massive scalar field. 

Grishchuk and Sidorov [26] discuss the initial conditions for the Wheeler-DeWitt equation, especially for the massive 
scalar field in a closed Friedmann model. 

In Ref. [27], Amendola, Khalatnikov, Litterio, and Occhionero consider quantum cosmology with a complex field. 

Guendelmann and Kaganovich [28] discuss cosmic time in quantum cosmology. The factor-ordering problem is 
solved such that the kinetic term gets the form □ + £R where £ < £ c <m/- 

A comparison of the minisuperspacc of minimally and conformally coupled scalar fields was done by Page [29]. He 
solves the factor-ordering problem of the Whcclcr-DcWitt equation by requiring that the kinetic term is proportional 
to the Laplacian in the minisuperspace metric, i.e., our variant B with e = 1. The classical equation is similar to the 
geodesic equation in superspace. 

Refs. [30] deal with quantum cosmology from the path integral point of view: The Wheeler-DeWitt equation can 
be derived in a first approximation from the corresponding path integral. Halliwcll [30] solves the factor-ordering 
problem by requiring invariance with respect to field redefinition of both 3-metric and lapse function. Jafariztadeh, 
Darabi and Rastegar [30] apply the method of Duru and Kleinert to evaluate the path integral for quantum cosmology, 
cf. Kleinert [30]. 

DeWitt [30] assumes the path integral to be the more fundamental approach, but the Wheeler-DeWitt equation in 
the minisuperspace to remain a good approximation to it. 

By considering the Wheeler-DeWitt equation 



Gibbons and Grishchuk [31] obtain the result that inflation is typical in the set of spatially flat Friedmann models 
in Einstein's theory with a A-term. The parameter p is due to the factor-ordering ambiguity, they take p = 1 as a 
preferred value. An analogous result is given by Hawking and Page [32] . Melnikov and Pevcov [33] discuss the factor- 
ordering problem for the Wheeler-DeWitt equation in closed and open Friedmann models and give some solutions to 
it. 

Reuter and Schmidt [34] and Schmidt [35] derive the Wheeler-DeWitt equation of fourth-order gravity for a spatially 
flat Friedmann model and compare with the corresponding conformally equivalent (due to the Bicknell theorem) second 
order models. The solution of Reuter/ Schmidt for flat Friedmann models is generalized in Ref. [36] by Pimentel and 
Obregon to closed and open models. Fabris and Reuter [36] continue to generalize the results of [34] to show that the 
Bicknell theorem applies also at the level of the Wheeler-DeWitt equation. 

Rainer [37] gives an overview on higher dimensions and discusses 3 types of conformal transformations of different 
levels for the Wheeler-DeWitt equation. Quite recently, the solutions of the Wheeler-DeWitt equation in comparison 
with the appearance of singularities is treated by Mongan [38] . Zhang and Shen [39] consider quantum cosmology 
with a complex scalar field at finite temperature. 
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A critical discussion of the Wheeler-DeWitt equation and an alternative quantization scheme is presented in Ref. [40] 
by Peres. 

From recent constraint calculations of Hwang and Noh [41] to an inflation model based on a non-minimally coupled 
massive scalar field and comparisons with observational data (COBE-DAIK), one can state that minimal coupling is 
a good approximation for inflationary models. Following Futamase and Maeda [42], the coupling constant is either 
quite small, £ < 1/1000, or negative. According to the coupling factor (1 — £0 2 ) in front of R, critical behaviour 
appears for £ = 0~ 2 , i.e. for positive values of £ only. 

Let us now summarize our results: We solved the Wheeler-DeWitt equation for the minisuperspace of a cosmological 
model of Bianchi type I with a minimally coupled massive scalar field 4> as source by generalizing the calculation of 
Lukash and Schmidt [1] . Contrarily to other approaches we allowed strong anisotropy. 

Combining analytical and numerical methods, we applied an adiabatic approximation for 0, and as new feature 
we found a period-doubling bifurcation of the typical solutions. This bifurcation takes place near the cosmologi- 
cal quantum boundary, i.e., the boundary of the quasiclassical region with oscillating ^-function where the WKB- 
approximation is good. The numerical calculations suggest that such a notion of a "cosmological quantum boundary" 
is well-defined, because sharply beyond that boundary, the WKB-approximation is no more applicable at all. 

This result confirms the adequateness of the introduction of a cosmological quantum boundary in quantum cos- 
mology. We applied the supercovariance principle, i.e., the underlying theory should also be covariant with respect 
to transformations representing a mixture between space-time and matter degrees of freedom. With our figures, we 
tried to vizualizc the birth of the Universe. 
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Until Eq. (4.6) we parallely dealt with variant A (e = 0) and B (e = 1), but beginning from Eq. (4.7) we simplified 
by restricting to variant A. Now we want to complete the calculation by showing what changes using variant B. 

The purpose of this appendix is to show in detailed calculations what has been verbally mentioned in the literature 
several times, namely the fact that differences in solving the factor-ordering problem do not essentially change the 
results. For this consideration we now set e = 1 in Eq. (4.6) and use the definitions below that equation: 
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APPENDIX: VARIANT A VS. VARIANT B 




(5.1) 



The transformation from Eq. (4.7) to Eq. (4.8) now leads to 




A 2 



)y = o, 



(5.2) 



whose solutions are also Bessel functions, modified by a factor x 



.-1/2. 



y{x) = x 




(5.3) 
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where C±, C2 are constants. 

This means that the wavefunction (4.17) using variant A is changed by multiplying by a factor ^/3he~ 3a /mip. 
Moreover, the index A of the Bessel functions must be replaced by -J A 2 + 1/4. 
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